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Abstract 



We study the direct CP violation induced by inelastic final state interaction (FSI) rescattering 
in _D — s- TTTT modes, and find that the resultant CP asymmetry is about 10~^ which is larger than 
e' in the K-system. Our estimation is based on well-established theories and experiment measured 
data, so there are almost no free parameters except the weak phase (5i3 in the CKM matrix. 
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I. Introduction. 



The study on CP violation phenomena is so important for understanding the mechanisms 
in particle physics that the subject has attracted great attention of both experimentalists and 
theoreticians of high energy physics for several decades already. The first and so far the only 
observation of CP violation is the measurement of e in the neutral K-system, as e ~ 2.3 x 10^'^ 
the non-zero e is due to a mixing of — K^ through box diagrams, which is theoretically evaluated 
at the quark level and the result predicts an approximate mass of charm-quark However, the 
parameter e' which is related to processes as a direct CP violation has not been reliably measured 
yet, it is believed that e'/e is non-zero and is expected between 10~* to 3 x 10"'^ 

In the well-known theory, the mechanism for resulting in non-zero e' is due to an interference 
of two pions in the neutral K-decays Q, and all the details are given in ref.Q. Since 27r can be 
in both 1=0 and 2 isospin eigenstates, if the two isospin channels have different weak and strong 
phases, their interference would result in a CP asymmetry proportional to sin((5o — 82) sin(0o ~ ^2) 
where 5i is the weak phase and Oj corresponds to the the phase shifts emerging in the strong 
interaction rescattering processes of 1=0 and 2 channels respectively. 

Generally, a direct CP violation must be realized through interferences between at least two 
amplitudes having different weak and strong phases, no matter at quark or hadron levels. The 
weak phase is determined by an underlying theory, for example, the CKM theory two-Higgs- 
doublet model |^ etc., while the strong phase can be either produced in the strong scattering at 
hadron level as in K-system or occur at the quark level. If an absorptive part of loops exists 
the strong phase is non-zero. For instance, at the high energy processes, B-physics [Q, top-physics 

and high energy collisions , the main part of the strong phases come from the absorptive part 
of loops, even though at the hadronization process, the hadron rescattering can also cause a strong 
phase. 

As the observation concerns hadron products, the strong final state interaction is not negligible 
and CP asymmetry may occur due to the phase shifts in the rescattering. Recently, by studying 
single pion exchange inelastic FSI for D — > VP processes, we pointed out that the inelastic strong 
FSI due to t-channel particle exchange may play important roles for producing CP violation in 
D and B hadronic decays A very recent estimation by Blok, Gronau and Rosncr shows 
that the inelastic FSI for B tttt, KK may produce th CP asymmetry as large as 10 ~ 20%, 
disregarding time dependent B-B mixing effects. Due to the GIM mechanism, the — mixing 
is small compared to B'^ — B'^ even though models beyond three-generation standard model 

may result in larger mixing effect | |15| . Recently, Browder and Pakvasa JlGf re-studied experimental 
implications of large CP violation and final state interactions (FSI) in a search for D° — _D" mixing 
and they concluded that FSI is important. Close and Lipkin studied possible strong FSI due 
to exotic resonances in D exclusive decays. They constrain their analysis at the most Cabibbo 
favored channels where the rescattering is elastic only. Obviously, it would be interesting to study 
CP asymmetry in D-system due to the inelastic strong FSI. 

In Z? — > TTTT decays, there are both elastic and inelastic rescatterings, in general, an amplitudes 
of D ^ TTTT should be 

T{D f) = r*''""(L' ^ /) + T^^\D /) (1) 

and 

T^si f\T\n> pn<n\H,ff\D> (2) 

n 

where |n > is a complete set of the strong interaction states satisfying four-momentum conservation, 
and pn is the density of state \n >. At ^/s = Md energy region, the most important intermediate 
states are TTTT, KK, pp, K*K*, aiir, a27r, KiK. These intermediate mesons can be on their mass 
shell, thus to be real particles. The off-shell contribution can be attributed into the quark level 
because the intermediate particles are virtual. 

In our present paper we only consider the elastic tttt — > tttt and the inelastic KK — > tttt 
rescattering processes which have experiment measurements [|l8| at the energy of D meson mass 
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and therefore are well constrained. This is different from the B meson case, where the inelastic 
scattering tttt ^ KK amplitudes can be only estimated by theoretical models . Contributions 
from other intermediate states may modify the results by adding a factor round unity, but are 
unlikely to change the whole scenario and order of magnitude of the CP asymmetry. 

Considering the direct tree level transition amplitude, tttt — s- tttt and KK tttt FSI, we have 

T{D tttt) EE ^ tttt) + T^'^^iD ^ tttt ^ tttt) + T^^\D KK tttt). (3) 

Here only the inelastic rescattering from KK intermediate states can induce a direct CP violation, 
but not the elastic ones. Since the tree amplitude of tttt and KK have different weak phases 

^"'"gi^cdVud) and Arg{V*gVus) respectively and the phase shifts in the inelastic rescattering KK ^ 
TTTT is non-zero, the interference between the two parts T*^'^'^{D tttt) and T^^'{D KK tttt) 
would result in a non-zero CP violation. The contribution of elastic scattering does not change 
the weak phase of T*^'^'^{D — > tttt), but they can cause a strong phase shift to T*'"'^'^. Namely, 

ftree^jj _^ ^ j^tree^j^ _^ ^ T^SI^j^ TTTT TTTt) = T^"""" X f 

where / is a scattering probability amplitude. 

Moreover, both D KK and Z? — > tttt are Cabibbo suppressed modes, so their tree amplitudes 
have the same order of magnitude. Even though, one expects that the FSI may change their relative 
ratios somehow JTol , the order remains the same. Their interference may be large, since two parts 
suffer the same Cabibbo suppression. 

Our numerical results show that the CP asymmetry can be about 10^"*. 

In next section, we present the formulation for evaluating the CP asymmetries, in the third 
section we give the numerical results and the final section is devoted to our conclusion and discus- 
sion. 



II. Formulation 

(i)The tree level amplitudes. 

The effective Hamiltonians for c ^ d + u -\- d and c ^ s + u -\- s are 



H 

and 



7/ = ^KrfKd[aid7M(l - 75)cS7'^(l - I5)d + a2dj^{l - j5)duj>'a - 75)0] + h.c, (4) 



H^fj = ^V:^Vus[aisj^{l - 75)cw7^(l - 75)s + 0257^1 - 75)s^l7^(l - 75)c] + h.c, (5) 
where the color indices are omitted and ai, 02 are parameters as 

Ol = Ci -f ^iC2, as = C2 + ^2Cl, (6) 

and ^ ^ 

where ri and r2 correspond to a non-factorizable contribution of the hadronic matrix elements 

ci = , and C2 = ^ ^ ' (8) 

where c± can be derived with the renormalization group equation, numerically at the energy scale 
of the charm quark, 

ci ~ 1.26, C2 ~ -0.51. 



2 



By fitting data, Cheng obtained ri ~ r2 ^ —0.67 for D PP decay Tlicn in om- later 

calculations, we use ai = 1.26 and a2 = —0.51. 

As many autliors suggested, we can ignore tlie contributions from the W-exchange and anni- 
hilation quark diagrams so the amplitude at the tree level can be obtained with the vacuum 
saturation approximation and the non-factorization effects are absorbed into the parameters ri 
and r2. We have 

<7r+n-\Hif^.\D°> = ^v:,VuaaM < 7r-\dj,{l - ^,)c\D" > 

= 9Zv;^Vu,a,Mf^^Ml-ml) + f^-ml], (9) 
where f^^ are the form factors in D to tt transition. With the multi-pole approximation p3| 



r2 



Uiq') = F,{q')=FM^)C7JH—^.T. (10) 



f-il') = (11) 



where Md* = 2.010 GeV is the nearest po le, qf^ == {Md — mT^Y . In the following, we take the 
single-pole approximation as n=l in eq.(|lO|). 

For — > TT^TT*^ and K^K^ we have similar expressions. Ignoring the W-exchange and 

annihilation, the tree amplitude for K^K^ is zero and the transition can only be realized 

through the elastic and inelastic FSI rescattering. 

The calculated /± ^ values have been checked by using the method given in ref. |2^ and found 
they coincide with each other very well. But Roberts' parameters [|4| are obtained by fitting the 
data oi D K transition only, so the results of /^'^ obtained in the multi-pole approximation 
deviate from that calculated in terms of the parameters of ref. |Q assuming an SU(3) symmetry. 
Therefore we later take only the values /^'^ obtained in the multi-pole approximation. 

(ii) The elastic and inelastic FSI. 

The S-matrix for strong interaction is 

where T-matrix is the non-trivial part determined by the strong interaction Lagrangian. 

It is noted that the 5mn term corresponds to a no-interaction scattering transition (or the 
trivial part of the S-matrix), so is exactly the "tree" part of eq.(j^). For the elastic and the inelastic 
rescattering contributions, the amplitudes read 

T^^^D"^ ^TT+TT- ^TT+TT-) = i < TT+TT" |T| TT+TT" > < TT+TT" liJ^^^^ l^" > , (13) 

T^'S/^^o ^0^0 ^ ^+^-) ^ z<^+^-|T|7rV>p, <7rO^O|ijg^^|7^o>, (14) 
T^^\D^ ^ K+K- ^TT+TT-) = i <Ti+TT-\T\K+K- > pK <K+K-\Hf^f\D^ > . (15) 
With help of the isospin analysis, for the elastic scattering, 

p. < ^+7r- ITK+TT- >= ^Toe'^" + ^T^e^'^ , (16) 

p. < 7r+7r-|T|^V" >= ^(-roe''^« +T2e*^^), (17) 

where Tq, 9q and T2, 62 are the measured scattering amplitudes and phase shifts of 1=0 and 1=2 
channels p8[, respectively. The transitions to the tt^tt" final state have similar expressions. 
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The contributions from the elastic FSI of tttt tttt can be absorbed into the tree amplitudes, 
in our case they do not provide a different weak phase from the tree amplitudes, but result in a 
strong phase shift. Including the elastic scattering, the amplitudes for D° tt+tt" and tt^tt" can 
be written as 

ftree^jjO ^ ^+^-) 



(18) 



(19) 

where the notation T*'"'^'^ refers to the tree amplitude modified by the elastic scattering and T,;, 6i 
are measured values. 

For the inelastic scattering, KK tttt, there are experimental measurements We can 
decompose K^K~ and tt+tt^ in terms of the basis of isospin as 

\K+K-> = ■^[\1,0>+\0,0>]k; (20) 



Thus 



k+^"> = [y-|o,o>+y-|2,o>],. (21) 

T" (1)° TT+TT-) = T^^^D"" K+K- n+TT-) 

^KlKs/xai [/f ^(Afl, - mi-) + f°^ml]z^TKe^''< ./^^ (22) 



where p^/p^ = ^1 - 4m2^/Af2 /^l _ Aml/M^. 

The similar expressions can be easily derived for T^^^{D'> ^ K+R- ^ tt^^O)^ 
(iii) The CP violation. 

As long as there are more than two channels with different weak and strong phases, their inter- 
ferences can result in a direct CP violation. If the total transition amplitude T is a superposition 
of two independent amplitudes Ai and A2 as 

T = Aie'^'e''!'' + ^26^*^6*"^% (23) 

while its CP conjugate amplitude is 

T = Aier'^'e"''' + Aae-'^^e*"^^ (24) 

Thus a direct CP asymmetry is defined as 



^ |r|2 + |r|2 

sin((5i - ^2) sin(0i - 02) 

(25) 



A( + Ai + 2A1A2 cos((5i - (52) cos(0i - (j)2) ' 
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In our case of D — > tt+tt (tt^tt"), the two interfering amplitudes are the modified "tree" part 
fitreep _^ Q ^j^g p^j.g "inglastic FSI" part r'"'^''^^"^ given in (|2|), which have different 
weak phases Si and strong phases For D'^ — s- tt'^tt'^ the expressions are similar. 
Let us work within the framework of CKM matrix 

VudV;^ - (ci2Cl3)(-Si2C23-Ci2S23Sl3e^''^^)* (26) 
VusV;, = (si2Cl3)(ci2C23-Si2S23Sl3e*^^^)*. (27) 



Thus 



6^-- = arctan -^i^^^3.i3 sin J13 . ^^g^ 

S12C23 + C12S23S13 C0Sdi3 

S^-^ = arctan ^i^^^3.i3 sin<5,3 ^ ^^g^ 

C12C23 — S12S23S13 COS(3i3 

In f*^<=^{D^ vr+TT-) = T^^'^'^iD" -> n+n-) + TP^^{D° tttt ^ tt+tt-), the weak phase 
is that of the tree amplitude as S"^"" , while that of T^^'{D° K+R- tt+tt-) is S^^^ , 
therefore they are different. One can also notice that |(5^^^| <C |(5^^'^| in this convention, even 
though the final result is independent of the convention adopted in the calculations. 

Later we will evaluate 

^ r(g"-.^+7r-)-r(5"->7r+^-) . 

^ r(LlO^ TT+TT-) +r(L)0^ TT+TT-)' ^ ^ 

^ r(D°^ A°)-r(^°-> A°) 
' " r(i?o ^ ^OttO) + r(^o ^0^0) ■ ^ ' 

In next section, we will present our numerical results of Ri and i?2- 



III. The numerical results 

(i) All the CKM entries involved in our calculations have been measured, even though there 
are some uncertainties 01 . Thus with the optimistic sin5i3 = 1, we have 

gD^^ ~ 3.2 X lO"'^ - 1.1 X 10"^ 6^^'^ ~ 5.3 x 10"^ 

Later we will use the most favorable values for the CP violation calculations. 

(ii) We obtain in terms of the dipole approximation, 

f^^q^^ ml) ^0.8; f^^q^= ml) ^-0.8; 

= O « 0.7; f^^iq' = ml) « -0.65. 

(iii) For the elastic and inelastic scattering tttt ^ tttt, KK ^ tttt, the transition probability 
amplitude T and the phase shift 9 are experimentally measured |jl^ as 

ro(TTTT ^ tttt) w -0.48, and 6*0 « 308°, 

T2(tttt ^ tttt) w -0.45, and 6*2 « -50°; 

T'"''^'"''"' = Tk{KK ^ tttt) « 0.1, and OKiKK tttt) « 310°, 

for the energy range of Md- So by the notation of eq.(p5|) sin((5i — ^2) ^ —1.1 x 10"^^ while 
02 ~ 310°, but 01 of eq.(p5|) must be evaluated by eg. (p^) and (p9[), in our case. 

(iv) The CP asymmetries. 

With the information given above, we obtain 

Ri = -1.1 X 10""^; (32) 
i?2 = 2.2 X 10"^ (33) 
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In these calculations we almost do not have any free parameters, except the CKM phase S13. We 
have taken sinSi^ = 1 and neglected the contribution from other intermediate states. They may 
modify the results, but are unlikely to change the whole scenario and order of magnitude of the 
CP asymmetry. It is noted that i?i and R2 have opposite signs, as a matter of fact, due to the 
uncertainty of 613 in the CKM matrix, the absolute sign of Ri is not important, but only the 
relative sign is meaningful. 

IV. Conclusion and discussion 

In this work, we discuss the CP violation effects caused by the inelastic FSI rescattering in 
D — > TT+TT^, tt'^tt'^ modes and obtain the CP asymmetry ratios of order 10~^. 

(i) The observed indirect CP violation in the K-system which is characterized by e is of order 
10""^. Even though the direct CP violation e'/e has not been reliably measured yet, present data 
incline to confirm it to be (1.5 ± 0.8) x 10"'^. Anyhow, one has all reasons to believe it is non- 
zero and of order 10^** ^ 3 x 10^'^ namely the superweak mechanism is almost ruled out by 
experiments. Our estimation on D —^ tttt, KK shows that the direct CP asymmetries here are 
about 2 orders of magnitude larger than e'. 

(ii) In the expression, one can see that the direct CP violation is caused by the interference 
between the tree amplitude modified by the elastic rescattering (denoted as T*^^^ in this paper) 
and that induced by the inelastic FSI rescattering while the two parts have different weak and 
strong phases. The interference is proportional to a product of the two amplitudes |'7"*''ee| jjiFS/j 
and the differences of weak and strong angles | sin((5^^ — S^^) ■ sin(</)i — (f>2)\- There are two factors 
which suppress the CP asymmetry values. 

The first one is that in the framework of CKM matrix, the weak phase is about order 10"'^ 
which is independent of convention. 

The second suppression factor comes from the measured strong phase difference | sin(0i —02)1 = 
0.3 and amplitude Tk{KK tttt) - 0.1. 

As discussed in ref. [l9j, the final state interaction effects can be described in a hadronic 
triangle diagram and the absorptive part of the loop gives rise to a strong phase. In our previous 
work, we only estimated the absorptive part and stressed that the inelastic FSI is important 
in many processes, so that cannot be ignored. The real part is hard to be properly evaluated 
because of the ultraviolet divergence and the obtained results would depend on the renormalization 
scheme. In this work, as suggested in literatures [ p5| , we only deal with the FSI, namely the 
intermediate hadrons are real particles, i.e. on their mass shell, the dispersive part of the loop 
is small compared to the tree amplitude, in fact, its contribution is effectively absorbed into the 
phenomenological parameters ai, 02 which may slightly deviate from the values derived in terms 
of the renormalization group equation, in our case (Note, it is not always true). Thus, we use only 
the data directly obtained from corresponding experiments, so can avoid any ambiguity caused by 
theoretical uncertainties. 

(iii) Since the proposed channels to be observed are Cabibbo-suppressed, the decay rates could 
be smaller than the Cabibbo favored channels by sin^ 6c roughly. 

(iv) The branching ratio oi D KK;7T7t is about 2 x 10^"^. The production cross section cr of 
D°D" is measured at BEPC [|6|, 

a{D°D°) = 11.63 ± 1.1 n6, at the BEPC energy. 

Taking the most optimistic values evaluated in the framework of C-K-M theory, the number of 
events for observing the Cabibbo-suppressed decay channels D KK; tttt would be 

TV = L X 2 X 10"'"^ X CT X r X /, (34) 

where L is the luminosity, r is the measuring time period and f is the observation efficiency. 
For the proposed charm-tau factory, L can reach 10'^'^ cm~^ sec~^ , so 

TV w 7.3 X 10^ X / X n, 
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where n is the number of necessary years. Since the CP asymmetry is ranged about 10"**, so to the 
reasonable statistical level for observation of CP violation, N at least must be 10^ ^ 10^, it would 
need a charm-tau factory with luminosity of 10'^'^ cm~^sec~^ to run for 15 years. Even though this 
number is not much encouraging, but as suggested by Browder and Pakvasa p6[ |, if there is new 
physics which can provide with a larger weak phase, the observation becomes very possible. Even 
with this small CP asymmetry, there is still possibility to make the measurement. 

Therefore, for measuring the direct CP violation which is one of the main interests in the field 
of high energy physics, a high luminosity r-charm factory would be extremely helpful. 
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